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Abstract 


This paper deals with the evolution equation of a curve obtained as the sharp 
interface limit of a non-linear system of two reaction-diffusion PDEs. This sys¬ 
tem was introduced as a phase-field model of (crawling) motion of eukaryotic 
cells on a substrate. The key issue is the evolution of the cell membrane (in¬ 
terface curve) which involves shape change and net motion. This issue can 
be addressed both qualitatively and quantitatively by studying the evolution 
equation of the sharp interface limit for this system. However, this equation is 
non-linear and non-local and existence of solutions presents a significant analyt¬ 
ical challenge. We establish existence of solutions for a wide class of initial data 
in the so-called subcritical regime. Existence is proved in a two step procedure. 
First, for smooth ( H 2 ) initial data we use a regularization technique. Second, 
we consider non-smooth initial data that are more relevant from the application 
point of view. Here, uniform estimates on the time when solution exists rely on 
a maximum principle type argument. We also explore the long time behavior of 
the model using both analytical and numerical tools. We prove the nonexistence 
of traveling wave solutions with nonzero velocity. Numerical experiments show 
that presence of non-linearity and asymmetry of the initial curve results in a net 
motion which distinguishes it from classical volume preserving curvature mo¬ 
tion. This is done by developing an algorithm for efficient numerical resolution 
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1. Introduction 

1.1. Phase field PDE model of cell motility 

This work is motivated by a 2D phase field model of crawling cell motility 
introduced in m- This model consists of a system of two PDEs for the phase 
field function and the orientation vector due to polymerization of actin filaments 
inside the cell. In addition it obeys a volume preservation constraint. In [T] this 
system was rewritten in the following simplified form suitable for asymptotic 
analysis so that all key features of its qualitative behavior are preserved. 

Let fl C K 2 be a smooth bounded domain. Then, consider the following phase 
field PDE model of cell motility, studied in Pj 

^=Ap e -±W'(p e )-P e -Vp e + \ e (t) in fl, (1.1) 

dP 1 

= eA P e - -P e - /3Vp s in fl, (1.2) 

where 

M*) = p- (jp 2 W \Ps) + P e ' dx, 

is a Lagrange multiplier term responsible for total volume preservation of p e , 
and 

W{p) = \p 2 {l~p) 2 (1.3) 

is the Allen-Cahn (scalar Ginzburg-Landau) double equal well potential, and 
(3 > 0 is a physical parameter (see [21]'). 

Q 



Figure 1: Sketch of phase-field parameter p £ 
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In (1.11-( 1.21, p e : Q — > R is the phase-field parameter that, roughly speaking, 
takes values 1 and 0 inside and outside, respectively, a subdomain D{t) C 
occupied by the moving cell. These regions are separated by a thin “interface 
layer” of width O(e) around the boundary T(t) := dD(t ), where p e (x,t ) sharply 
transitions from 1 to 0. The vector valued function P e : 11 — > R 2 models the 
orientation vector due to polymerization of actin filaments inside the cell. On 
the boundary 911 Neumann and Diriclilet boundary conditions respectively are 
imposed: d v p e = 0 and P e = 0. 

It was shown in jl] that p e (x,t) converges to a characteristic function XD(t) 
as e —»• 0, where D(t) C R 2 . Namely, the phase-field parameter 0 < p E {x, t) < 1 
is equal to 1 when x € D{t) and equal to 0 outside of D{t). This is referred to 
as the sharp interface limit of p e and we write T(t) := dD(t). More precisely, 
given a closed non self-intersecting curve T(0) C R 2 , consider the initial profile 

P ,(x.o ) = « c ( dist(l ; r(0)) ) 


where 


»„(,) = i (tank 


+ 1 


is the standing wave solution of the Allen-Cahn equation, and dist(x, T(0)) is 
the signed distance from the point x to the curve T(0). Then, p e (x,t) has the 
asymptotic form 

fe( ^Q=o.( dist(l ; r(t)) )+ow. 

where T(t) is a curve which describes the boundary of the cell. 


It was formally shown in jl that the curves T(t) obey the evolution equation 

V(s,t) = K(s,t) + ^(V(s,t))--^^- J K(s',t) + ${V(s',t))ds', (1.4) 

where s is the arc length parametrization of the curve T(t), V(s, t) is the normal 
velocity of curve T(f) w.r.t. inward normal at location s, |T(i)| is the length of 
T(t), k(s,1) is the signed curvature of T(i) at location s, and $(•) is a known 
smooth, non-linear function. 


Remark 1 . In JTj it was shown that 3>(V) = /3 ( E , o(I A ) where $o(V) is given by 
the equation 

MV) ■■= [ ^z;V)(6' 0 ) 2 dz (1.5) 

JR 

and 'ijj(z) = V) is the unique solution of 


^"(z) + Vtp'iz) - i/j(z) + 9’ 0 = 0, 


( 1 . 6 ) 


with ^>(±oo) = 0. 
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The case [3 = 0 in equations leads to = 0 in ( |1.4[ ), thus reducing 

to a mass preserving analogue of the Allen-Cahn equation. Properties of this 
equation were studied in sun], and it was shown that the sharp interface limit 
as e —> 0 recovers volume preserving mean curvature motion: V = k— A f r nds. 


Equations (1.1)-(1.2) are a singularly perturbed parabolic PDE system in 


two spatial dimensions. Its sharp interface limit given by (1.4) describes evo¬ 
lution of the curve T(t) (the sharp interface). Since V(s,t) and n(s,t) are 


expressed via first and second derivatives of r(f)(= T(s,t)), equation ( |1.4| ) can 
be viewed as the second order PDE for P(s,f). Since this PDE has spatial di¬ 
mension one and it does not contain a singular perturbation, qualitative and 
numerical analysis of (1.4) is much simpler than that of the system <0>-& 

Remark 2 . It was observed in m that both the analysis and the behavior of 
solutions of system (1.1)-(1.2) crucially depends on the parameter (3. Specifically 
there is critical value j3 cr such that for (3 > (3 cr complicated phenomena of non¬ 
uniqueness and hysteresis arise. This critical value is defined as the maximum 
/3 for which V — /34>o(E) is a monotone function of V. 

While this supercritical regime is a subject of the ongoing investigation, in 
this work focus on providing a rigorous analysis of subcritical regime (3 < f3 cr . 


For equation (1.4) the latter regime corresponds to the case of monotone function 
V-$(V). 

1.2. Biological background: cell motility problem 

In |2T] a phase field model that describes crawling motion of keratocyte cells 
on substrates was introduced. Keratocyte cells are typically harvested from 
scales of fish (e.g., cichlids [13]) for in vitro experiments. Additionally, humans 
have keratocyte cells in their corneas. These cells are crucial during wound 
healing, e.g., after corrective laser surgery PH- 

The biological mechanisms which give rise to keratocyte cell motion are 
complicated and they are an ongoing source of research. Assuming that a di¬ 
rectional preference has been established, a keratocyte cell has the ability to 
maintain self-propagating cell motion via internal forces generated by a protein 
actin. Actin monomers are polarized in such a way that several molecules may 
join and form filaments. These actin filaments form a dense and rigid network 
at the leading edge of the cell within the cytoskeleton, known as the lamellipod. 
The lamellipod forms adhesions to the substrate and by a mechanism known 
as actin tread milling the cell protrudes via formation of new actin filaments at 
the leading edge. 

We may now explain the heuristic idea behind the model. Roughly speaking, 
cell motility is determined by two (competing) mechanisms: surface tension and 
protrusion due to actin polymerization. The domain where p e {x) ss 1 is occupied 
the cell and P e as the local averaged orientation of the filament network. Surface 


double-well potential (1.3): 


tension enters the model (1.1)-(1.2) via the celebrated Allen-Cahn equation with 




4 























In the sharp interface limit (e —> 0), surface tension leads to the curvature driven 
motion of the interface. The actin polymerization enters the system (1.1)-(1.2) 
through the — P e ■ V/? e term. Indeed, recall 


Dpe 

Dt 


dpe 

dt 


Pe ' Vp £ 


as the material derivative of p E subject to the velocity field P e . Thus the term 
— P e ■ V p e is an advective term generated by actin polymerization. 


The last term of (1.1), A e (t) is responsible for volume preservation, which is an 


important physical condition. The diffusion term eAP e corresponds to diffusion 
of actin and does not significantly affect the dynamics of p e . The term —/3Vp E 
describes the creation of actin by polymerization, which leads to a protrusion 
force. It gives the rate of growth of polymerization of actin: ~ —/3Vp e . 

The Ap e term provides decay of P £ away from the interface, for example due to 
depolymerization. 

The system ( |1.2[ ) is a slightly modified form of the model proposed in 
m ■ It preserves key features of the qualitative behavior yet is more convenient 
for mathematical analysis. 

1.3. Overview of results and techniques 

A main goal of this work is to prove existence of a family of curves which 


evolve according to the equation (1.4) (that describes evolution of the sharp 


interface) and investigate their properties. The problem of mean curvature type 
motion was extensively studied by mathematicians from both PDE and geom¬ 
etry communities for several decades. A review of results on unconstrained 
motion by mean curvature can be found HEnna. Furthermore the viscos¬ 
ity solutions techniques have been efficiently applied in the PDE analysis of 
such problems. These techniques do not apply to mean curvature motion with 
volume preservation constraints EUSlEj, and the analysis becomes especially 
difficult in dimensions greater than two 0. Note that existence in two dimen¬ 
sional mean curvature type motions were recently studied (e.g., dm El) and 
appropriate techniques of regularization were developed. Recently, analogous 
issues resurfaced in the novel context of biological cell motility problems, after 
a phase-field model was introduced in |21| . 

The problem studied in the present work is two dimensional (motion of 
a curve on the plane). The distinguished feature of this problem is that the 
velocity enters the evolution equation implicitly via a non-linear function V — 
<I>(TA). Therefore the time derivative in the corresponding PDE that describes 
the signed distance, u(a, t) of the curve from a given reference curve also enters 
the PDE implicitly, which leads to challenges in establishing existence. 

The following outlines the basic steps of the analysis. First, we consider 
smooth ( H 2 ) initial data and generalize the regularization idea from [5] (see also 
0) for the implicit dependence described above. Here, the main difficulty is to 
establish existence on a time interval that does not depend on the regularization 
parameter e. To this end, we derive L 2 (in time and space) a priori estimates 
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independent of e for third order derivatives and uniform in time L 2 estimates 
for second order derivatives. These estimates allow us to show existence and to 
pass to the limit as e —> 0, and they are derived by considering the equation 
for u a = We use “nice” properties of this equation to obtain higher order 
estimates independent of e (which are not readily available for the equation 
for u). In particular, it turns out that the equation for u a can be written 
as a quasi linear parabolic PDE in divergence form. For such equations quite 
remarkable classical results establish Holder continuity of solutions for even for 
discontinuous initial data f 141 . This provides a lower bound on the possible blow 
up time, which does depend on H 2 norm of initial data for u in our problem. 

As a result, we establish existence on a time interval that depends on the 
H 2 norm of initial data. 

Second, observe that experiments for cell motility show that the cell shape is 
not necessarily smooth. Therefore one needs to consider more realistic models 
where smoothness assumptions on initial conditions are relaxed. In particular, 
one should allow for corner-type singularities. To this end, we pass to generic 
W 1,oc initial curves. For the limiting equations for u and u ff , we show existence 
on a time interval that does not depend on the H 2 (and H 1 for u a ) norms of 
initial conditions. This is necessary because these norms blow up for non-smooth 
functions from W 1 ’ 00 \ H 2 . The existence is proved by a maximum principle 
type argument, which is not available for the regularized equations that contain 
fourth order derivatives. Also it is crucial to establish Holder continuity results 
for u a , rewriting the equation for u a as a quasilinear divergence form parabolic 
PDE. 


After proving short time existence we address the issues of global existence 
of such curves. The latter is important for the comparison of theoretical and 
experimental results on cell motility. We will present an exploratory study which 
combines analytical and numerical results for the long time behavior of the cell 
motility model. Analytically, we prove that similarly to the classical curvature 
driven motion with volume preservation, traveling waves with nonzero velocity 
do not exist. While through numerical experiments, we observe a nontrivial 
(transient) net motion resulting from the non-linearity and asymmetry of the 
initial shape. This observation shows an essential difference from the classical 
area preserving curvature driven motion. 

Numerically solving (1.4) is a nontrivial task due to the non-linearity and 
non-locality in the formulation of the normal velocity. Classical methods such 
as level-set methods cannot be readily used here. We introduce an efficient 
algorithm which separates the difficulties of non-linearity and non-locality and 
resolves them independently through an iterative scheme. The accuracy of the 
algorithm is validated by numerical convergence. Our numerical experiments 
show that both non-linearity and asymmetry of the initial shape (in the form of 
non-convexity) can introduce a drift of the center of mass of the cell. Increasing 
the effects of the non-linearity or increasing the asymmetry results in an increase 
in the drift distance. 
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2. Existence of solutions to the evolution equation ( |1.4[ ) 

We study curves propagating via the evolution equation ( |1.4[ ). The case 
$ = 0 corresponds to well-studied volume preserving curvature motion (see, 
e.g., HHHHE]). We emphasize that the presence of <I>(E) results in an implicit, 
non-linear and non-local equation for V, which leads to challenges from an 
analytical and numerical standpoint. 

The goal of this section is to prove the following: 

Theorem 1. Let $ £ L°°(]R) be a Lipschitz function satisfying 


l$'l 


< 1 . 


( 2 . 1 ) 


Then, given T 0 £ W 1,0 °, a closed and non self-intersecting curve on R 2 , there is 
a time T = T(Tq) >0 such that a family of curves T(f) £ H 2 exists for t £ (0, T] 


which satisfies the evolution equation (1.4) with initial condition T(0) = Pq. 


Remark 3. The classes W 1,0 ° and H 2 above refer to curves which are parametrized 
by mappings from Sobolev spaces IT 1 ’ 00 and H 2 correspondingly. 

Remark 4. After time T the curve could self-intersect or blow-up in the parametriza- 
tion map (e.g., a cusp) could occur. 

We first prove the existence for smooth ( H 2 ) initial data. The main effort 
is to pass to non-smooth initial conditions (e.g., initial curves with corners), see 
discussion in Section [131 

Proof of Theorem [T] The proof of Theorem [T] is split into 4 steps. In Step 1 
we present a PDE formulation of the evolution problem (1.4) and introduce its 


regularization by adding a higher order term with the small parameter e > 0. 
In Step 2 we prove a local in time existence result for the regularized problem. 
In Step 3 we establish a uniform time interval of existence for solutions of the 
regularized problem via a priori estimates. These estimates allow us to pass 
to the limit e —> 0, which leads to existence for (1.4) for smooth initial data. 
Finally, Step 4 is devoted to the transition from R 2 -smooth initial data to W 1,00 
ones. A crucial role here plays derivation of L°° bounds for the solution and its 
first derivative independent of H 2 norm of the initial data. 


Step 1. Parametrization and PDE forms of (1.4). 


Let fo be a C 4 smooth reference curve in a small neighborhood of To and let 
To be parametrized by arc length parameter a £ I. Let fto(cr) be the signed cur¬ 
vature of To and v(a ) be the the inward pointing normal vector to To. Consider 
the tubular neighborhood 

Us 0 ■■= {x £ R 2 | dist(:r,fo) < 26 0 }- 

One can choose Tq and Sq so that the map 


Y: I x (—2(5 0 , 2<5 0 ) ->• U So , Y(a,u) := fo(cr) + uiz(a) 
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is a C 2 diffeomorphism between I x (— 2do,2do) and its image, and To C Y(I x 
(— 2Sq, 2do)) ■ Then To can be parametrized by To = To((t) + uo(cr)iy(a), a £ I, 
for some periodic function iio(a). Finally, we can assume that do is sufficiently 
small so that 

5o||koIU» < 1, (2-2) 

where kq denotes the curvature of F. 

A continuous function u : I x [0,T] —> [—do,do], periodic in the a variable, 
describes a family of closed curves via the mapping 

r(cr, t) = f 0 (cr) + u(a, t)v(a). (2.3) 


That is, there is a well-defined correspondence between T(cr, t) and u(cr,t). 
Recall the Frenet-Serre formulas applied to JV 


dr 

da 


(<r) 


n 0 (a)v{a) 


(2.4) 


where r is the unit 
velocity V of V as 


^(cr) = -K 0 (cr)r(cr) 
ent vector. Using (2.3)-(2.5) we express the 

1 — UKq 


(2.5) 
normal 

( 2 . 6 ) 


where S = S(u) = \/ii* + (1 — un o) 2 . Also, in terms of u, curvature of T is 
given by 

k(u) = Jg ((1 - UKo)u ac r + 2 K 0 ul + {K 0 ) a U a U + K 0 (l - UK 0 ) 2 ^ . (2.7) 


Note that if u is sufficiently smooth [4] one has 







in particular this holds for every u £ H per {I) such that \u\ < So on I. Here¬ 
after, H per {I) denote the Sobolev spaces of periodic functions on / with square- 
integrable derivatives up to the k-th order and W p £°(I) denotes the space of 
periodic functions with the first derivative in L°°(I). 


Combining (2.6) and (2.7), we rewrite (1.4) as the following PDE for u: 


u t - 


1 — UK 0 




1 — UKo 


Ut = 



Sda 


2t rj 
( 2 . 8 ) 


where 


L[u\ = J S{u)do 


is the total length of the curve. 

The initial condition u(a , 0) = uo(cr) corresponds to the initial profile T(0) = 
T 0 . Since (2.8) is not resolved with respect to the time derivative Ut , it is natural 


to resolve (1.41 with respect to V to convert (2.8) into a parabolic type PDE 
where the time derivative u t is expressed as a function of u, u a , u aa . The 
following lemma shows how to rewrite (2.8) in such a form. This is done by 


resolving (1.4) with respect to V to get a Lipschitz continuous resolving map, 


provided that £ L° 


and 


< 1 . 


It is convenient to rewrite in the form (1.41 


V = k{u) + 4 > (E) — A, J VS{u)d<r = 0, 


(2.9) 


where both normal velocity V and constant A are considered as unknowns. 

Lemma 1. Suppose that $ G L°°(R) and ||4> , || i oo( R ) < 1. Then for any 
Tt(cr) G Hp er (I), such that |zt(cr)| < Jq, there exists a unique solution (V(a), A) G 


L 2 (J) x R of ( |2.9[ ). Moreover, the resolving map T assigning to a given u G 
Hper(I) ^ { u i M — M the solution V = T(u) £ L 2 (I) of (2.9) is locally Lips¬ 
chitz continuous. 


Proof. Let J := ||$'|| L ° 


Fix k G C per (I), a positive function S G C per (I) 


and AgK, and consider the equation 


V = k + <f>(V) - A. 


( 2 . 10 ) 


It is immediate that the unique solution of ( |2.10 1 is given by V = ^(k — A), 
where is the inverse map to V — $(14). Note that yA_ < T' < therefore 
T is strictly increasing function and 'P(k; — A) —► ±oo uniformly as A —> =Foo. It 
follows that there exists a unique A G R such that V = 4/(k — A) is a solution 


of (2.10) satisfying 


VSdcr = 0. 


( 2 . 11 ) 
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V' 2 Sd(j 


1/2 


Next we establish the Lipschitz continuity of the resolving map (k, S) K > V £ 
L 2 (I) as a function of k £ L 2 (I) and S £ L™ r (I), S > 1 — <5o||ko||l°°(j) > 0 (cf. 
(2.2)), still assuming that k, S £ C per (I). Multiply (2.10) by VS and integrate 
over / to find with the help of the Cauchy-Schwarz inequality, 

J V 2 Sda = J nVSda + <S>(V)VSdo 

< I|5 , ||l«(7)(||k|U*(7) + (|/|||$||7~(K) 

Recalling that S > w := 1 — 5o||ko||l°°( 7) > 0, we then obtain 

II^11 t, 2 (7) < (||«||z, 2 (7) + (1^1 II < J > IIt j °°(m)) 1/2 ) - ( 2 -12) 

To see that k ^ V (for fixed S) is Lipschitz continuous, consider solutions Vi, 
Vi of (2.10)-(2.11) with k = «i and k = k 2 . Subtract the equation for V 2 from 
that for V[ multiply by S (Vi — V 2 ) and integrate over I, 

J((V 1 - V 2 ) 2 - ($(Vi) - $(V 2 ))(Vi - V 2 ))Sda = J( Kl - k 2 )(Vi - V 2 )Sda. 

Then, since |$(Vl) — $(V 2 )| < J|Vi — V 2 \ we derive that 

||^L - V 2 11 Zy 2 (7) < Ji/2(1 -J)W Kl ~ K ^ L2 V)- 


Next consider solutions of (2.10)-(2.11), still denoted V\ and V 2 , which corre¬ 
spond now S = Si and S = S 2 with the same k. Subtract the equation for V 2 
from that for V± multiply by Si (Vi — V 2 ) + V 2 (Si — S 2 ) and integrate over I to 
find 

J{{Vi - V 2 ) 2 -($(Vr) - $(V 2 ))(V! - V 2 j)S 1 da 

= J V 2 {S\ - S 2 )($(Vi) - $(V 2 ) - (Vi - V 2 ))da. 


Then applying the Cauchy-Schwarz inequality we derive 


1 + J 


||Vi - V 2 \\ L 2(!) < — — — 1|V 2 || L 2 (J) ||S : 


-s, 


2II 7°° (7) ■ 


Thanks to (2.12) this completes the proof of local Lipschitz continuity of the 
resolving map on the dense subset (of continuous functions) in 

e = {(K, S) £ L 2 (I) x {S £ L“ r (J); S > a,}} 

and thus on the whole set 0. It remains to note that the map u <—> (n(u), S(u)) 
is locally Lipschitz on {u £ H 2 (I); |u| < So on /}, which completes the proof 


of the Lemma. 


□ 
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Remark 5. The parameter A £ R with the property that the solution V of 
(2.101 satisfies f I VSda = 0 is easily seen to be 


A = 


L[u 


j (k(u) + $(V)Sd<j = 


7 r 

L[u\ 


1 

L[u] 


&(V)Sdcr. 


(2.13) 


Remark 6. Under conditions of Lemma 0 if k{u) £ H^ er {I) (u £ Hp er (I)) 
then it holds that V = F(u) £ Hp er (I). 


as 


Equation (1.41 is equivalently rewritten in terms of the resolving operator T 

S{u) 


u t = 


1 — UKq 


T(u), or ut = 


(2.14) 


where T(u) := S(u)J 7 (u)/( 1 — ukq). 

Step 2. Introduction and analysis of regularized PDE. 

We now introduce a small parameter regularization term to (2.8) which 
allows us to apply standard existence results. To this end, let u e = u E (a,t ) 
solve the following regularization of equation (2.14) for 0 < e < 1, 


+ eu laaa = ^(u E ), 


(2.15) 


with u E (a, 0) = Uq. Define 


rrF 1 — U E KQ . p , , 

V '■= -^- K + «W<r) ) 

where S e = S(u e ). Since V e = F{u E ), then by definition of the resolving map 
T we have that u e satisfies the following equation: 


■£U n 


ce q i 

$(W) = 


1 — UK o 


1 — U £ Kq 

S e 


k(u e ) 


(1 — U e Kq)L[u e ] 


^J^(V e )S E da + 2ny 


(2.16) 


Hereafter we consider H 2 er (I) equipped with the norm 


I H 2 (I) 


<L 2 (I) 


I L 2 (I)' 


Proposition 1. Let To and 4>(-) satisfy the conditions of Theorem^ 7] Assume 
that u o £ Hp er and max |ito| < Jo- Then there exists a non-empty interval [0, T E ] 
such that a solution u E of (2.16) with initial data u E (r],0) = uo(rf) exists and 

u e GL 2 (0,T e - H* er (I)) D R 1 (0, T E ; L 2 (I)) D L°°(0, T E \H 2 er ) 


and sup ||-u e (t)|| i oo (/) < J 0 - 
te[o,T] 


(2.17) 


Furthermore, this solution can be extended on a bigger time interval [0, T E + At] 
so long as 6 o — maxu(a,T e ) > a and ||m(T £ )||# 2 (/) < M for some a > 0 and 
M < oo, where A t depend on e, a and M, A t := A t (o, M, e) > 0. 
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Proof. Choose T £ > 0 and M > ||wo||/f 2 (/) ; and introduce the set 

K := {u : ||u||i /2 < M, |it| < <5o in I}. 

Given u £ L°°(0, T e ; K), consider the following auxiliary problem 

Ut P 7( j — F [vf (2.18) 

with u(cr, 0) = uo(a). Classical results (e.g., p”5]) yield existence of a unique 
solution u of (|2. 181 which possesses the following regularity 


u £ I/ 2 (0, T e ; Hp er (I)) fl ff 1 (0, T E ; L 2 {I)) fl L°°(0, T e ; H 2 er (I)). 

That is, a resolving operator 

T: L°°(0,T E -K) -> L°°(0, T e ; H 2 (I)) 

which maps u to the solution u is well defined. Next we show that F is a 
contraction in I \, provided that T E is chosen sufficiently small. 

Consider U\,U 2 £ T°°(0,T e ;/\) satisfying the initial condition iii(cx,0) = 
■U 2 (<U 0) = uo(cr) and define u\ := T(ui), 112 '■= T{u 2 )- Let u := u± — 112 . Then 
multiply the equality it t + £ii cr(T0 . cr = F(ui) — F(v, 2 ) by ( u ai7a(7 + u ) and integrate. 
Integrating by parts and using the Cauchy-Schwarz inequality yields 


1 d 

2 dt . 


^ (.^<7a u 2 )da + e || u <T<ra<T\\%2(l) + £ II u ctct||l2(7) 


'U'acracr II L 2 (I) + |M|l 2 (J))- 

Note that by Lemma [l] the map F(u) with values in L 2 (I) is Lipschitz on K- 
since F{u) = S(u)F(u)/(l — ukq) it is not hard to see that F(u) is also Lipschitz. 
Using this and applying Young’s inequality to the right hand side we obtain 


2 ~^.W u Wh 2 (i) - L'll'Wi - W2||^-2(j) + ^\\ u \\h 2 (i) 


(2.19) 


with a constant C independent of u\ and U 2 and T e . Applying the Gronwall 
inequality on (|2.19) we get 


sup ||T(iti) -T(u 2 )\\ 2 H 2 {i) < 2(e T - l)C||«i - u 2 \\l°°{p,T‘-,h*(i)) (2.20) 

0<t<T £ 


Similar arguments additionally yield the following bound for u = T(u), 


Sup |Hf)||^ 2( n 
0 <t<T £ 


< (e T - l)Ci + e 1 


W)’ 


( 2 . 21 ) 


with Ci independent of u £ L°°(0, T e ; K). Choosing T E sufficiently small we 
get that \\u(t) \\h 2 (i) < M for 0 < t < T E . 

Finally, multiply (2.181 by (it — Uq) and integrate. After integrating by parts 


and using the fact that ||it||jj 2 (j) < M we obtain 


sup \\u{t) - Mo||| 2 ( n < C 3 (e T - 1). 

0<t<T £ 
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Then using the interpolation inequality 


' — u o\\c(i) ^ C\\ u — uq\\h 2 (i)\\u — ito||i2(/) 


we get the bound 


sup \\u(t) - wollcm < C 4 II u(t) - «o||i 2( n < C 5 (e T ~ - 1). 

0<t<T £ 


( 2 . 22 ) 


Now by (2.20) and (2.22) we see that, possibly after passing to a smaller T E , 
T maps K into K and it is a contraction on I\. □ 

Step 3. Regularized equation: a priori estimates, existence on time interval 
independent of e, and limit as e —> 0. 

In this step we derive a priori estimates which imply existence of a solution of 
(2.16) on a time interval independent of e. These estimates are also used to 
pass to the e —> 0 limit. 

Lemma 2. Assume that uq € H 2 er and ||uo||l°°(j) < <$o- Let u e solve ( 2.16| ) on 
a time interval [0,T e ] with initial data u E ( 0) = u 0 , and let u E satisfy \u E (a,t)\ < 
Sq on I x T e . Then 


Il 2 (/) 


< a(t), 


where aft ) is the solution of 


d = 2Pa 3 + 2Q, a(0) = ||(VW|ll 2 


(2.23) 


(2.24) 


(continued by +00 after the blow up time), and 0 < P < 00 , 0 < Q < 00 are 
independent of e and uq . 

Proof. For brevity we adopt the notation u := u e , V :=V e , S := S E . Differen¬ 
tiate equation (2.16) in a to find that 

d 5 u 


Ufjt + £ 


da 5 


S &(V)V a = , S (k(u))„ 


1 — UKo 

s 


1 — UKo 


+ (-—) (*(«)-$(V) 

V 1 — UKo ' ° '■ 


1 


27T 


L[u\ L[u\ Jj 


$(V)Sda 


(2.25) 


Next we rewrite (2.16) in the form V = $(F) + k — A to calculate V a , 

(1 - &(V))V a = K a 

whence 


V a = 


1 


1 - $'(V) 

Now we substitute this in (2.25) to find that 

d U Uaaa 


Uat + £ 


$'(V)+2 

dcV “ 5 2 (1-$'(V)) ~ (1 - $'(V))S 4 

S 


UaU 2 aa + A(<J, V, U, U a )Uacr 


+ B(a,V,u,Ua) + (. -) ($00 

VI — UHlr) / <T V 


27r 


L[u] L[u\ 


$>{V)Sda), 

(2.26) 
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where A(a. V. u. p) and B(a,V,u,p ) are bounded continuous functions on / x 
R x [—(Jq, d 0 ] x R. 

Multiply (2.26) by u aar7 and integrate over /, integrating by parts on the 
left hand side. We find after rearranging terms and setting 7 := sup(l — $'(V)) 
(0 < 7 < 00 ), 


2 l 2 (/) + £ \\ u <™m\\ 2 L 2 (i) + - J | 


?%<C 


+ CA ||«„ 


I L 2 (I) 




| | ^aaa 


U, 


. da 

2 da\ 1/2 

aaa \ g 2 J > 

(2.27) 

where we have also used the inequality ||u<t|||oo(/) < 1 1 | ||iw||£ 2 (j) and esti¬ 
mated various terms with the help of the Caucliy-Schwarz inequality. Next we 
estimate the first term in the right hand side of (2.27) as follows 

, , 9 , , da 


c(7 1 | u 'CT(T(j | ^3 _ J | “'era \ 


l - + — f\ 

s 4 4 C 7 Jj 


, da 
52 


(2.28) 


and apply the following interpolation type inequality, whose proof is given in 
Lemma[3j for all u £ Hp er (I) such that |u| < do it holds that 

J\ u ^<»j\u^ + C 2 (Jul„do) + C 3 V/z > 0, (2.29) 

where C 2 and C 3 depend only on p. Now we use (2.28) and (2.29) in (2.27), and 
estimate the last term of (2.27) with the help of Young’s inequality to derive 
that 


1 d 

2 Jt" u ' 


1 f da 

wlli 2 (J) + £ \\ u a<T<Ta\\L 2 (I) + ^ J \u aiJr ,f-^ < P\\u aa \\ G L i^ + Q. (2.30) 


Then by a standard comparison argument for ODEs ||zw||| 2 ( 7 ) < a(t), where 
a solves (2.24). The Lemma is proved. □ 


Lemma 3. Assume thatu £ Hp er (I) and |it| < do on I. Then (2.29) holds with 
C 2 and C 3 depending on p only. 

Proof. The following straightforward bounds will be used throughout the proof, 


Note that 


7 ^(K| + (1 - do||Ko|| L »( 7 ))) < S{u) < \u a \ + 1, 


\(S(u)U < — {\UaWUaal + C\u a \ + C{). 

^ C \\ u la/S\\ Lao{I) Jul a da, 


(2.31) 

(2.32) 

(2.33) 
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where C > 0 is independent of u. Since f r u aa da = 0 we have 

\\u 2 a JS\\ L oo (I) < J^ulJSUda. (2.34) 


Next we use (2.31), (2.32) to estimate the right hand side of (2.34) with the 

da 


help of the Cauchy-Schwarz inequality, 

da 


J\{ul a /S) a \da < 2 J\u, 


Uac r I ry 4“ C 


( 7 ( 7(7 (^( 7(7 


< 2 


S 

1/2 


s 2 


il„da 


U La^j ( / 


1/2 


(2.35) 


+ C[ iuU%) 


1/2 


f^lado^ +C J^ul a da. 


Plugging (2.341-( 2.35) into (2.33) and using Young’s inequality yields 


da 


r *C^I<2 


da\ 


1/2 


u 2 a<J da 


+ 


s 2 J / 

da 


\ 3/2 


cUu^da^j + C. 


Finally, using here Young’s inequality once more time we deduce (2.29). 


D 


Consider a time interval [0, T*] with T* > 0 slightly smaller than the blow 
up time T bu of (2.24). As a byproduct of Lemma [$] (cf. ( |2.30| ) ) we then have 
for any 0 < T < min{T £ ,T*} 


pT pT 

sup ||W(7cr||i2/j^ + £ / ||lt(Tcro-(T||i 2 (/)dt + 

te[o,T] 'Jo 


J j haoaWll (I) dt<C (2.36) 


where u := u e (a,t) is the solution of (2.16) described in Proposi tion |T| and C is 
independent of e and T. In order to show that the solution of (2.16) exists on 


a time interval independent of e, it remains to obtain a uniform in e estimate 
on the growth of ||u £ — Mojjcm i n time. Arguing as in the end of the proof of 
Proposition [l] and using (2.36) one can prove 


Lemma 4. Assume that uq £ H 2 er (I), ||ito||z,”(/) < ® n d assume that the 


solution u e of (2.16) satisfies \u e (a,t)\ < Sq on I x [0,T £ ]. Then for all 0 < t < 
min {T £ ,T*}, 

Wuf-Uofc^^C^-l) (2.37) 

where C is independent of e. In particular, there exists 0 < T** < T*, indepen¬ 
dent of e, such that sup{ || u e (t )||(/)! 0<t<min{T £ ,T**}} < 5 0 . 

Combining Proposition [T] with Lemma [2] and Lemma [4] we see that the so¬ 
lution u e of (2.161 exists on the time interval [0,T**] and ( |2.36 ) holds with 


15 





















_ rj-y* 


Now, it is not hard to pass to the limit e 


0 in (2.16). Indeed, 
u weak 


Combining 


exploiting (2.36) we see that, up to extracting a subsequence, u‘ 
star in L°°(0, T**; H 2 er (I)) as e —> 0. Using (2.36) in (2.15) we also con¬ 
clude that the family {Mf}o<e<i is bounded in LF(I x [0, T" 
uniform estimates on ||A 
deduce u e 

itf +eu% aacr —*■ ut weak star in L°°(0, T**; L 2 {I)) and ztf + eu% aaty — > u t strongly 
in L 2 (0,T 


"gcjg IIl 2 (/x[o.t]) (from ( |2.36D ) and |K||l2(/ x [o,t**]) we 
u strongly in L 2 (0,T**; H 2 er (I)) D (7(0, T**; H) er (I)). Moreover, 

n** . r 2 / 

M GGGG fTV/MUA WVMU. ^ ^ -*■ 5 

n **; Hp er (I)). Thus, in the limit we obtain a solution u of (2.8). That 
is we have proved 


Theorem 2. (Existence for smooth initial data) For any uq € H 2 (I) such 


that ||wo||l“(/) < <5o there exists a solution u of (2.8) on a time interval [0,T], 
with T > 0 depending on uq. 

Remark 7 . Note that the time interval in Theorem^can by chosen universally 
for all uq such that ||itoHz, 00 (zj < a < <5o and ||itoIIzr 2 (z) < M < oo, that is 
T = T(M, a) > 0. 

Step 4. Passing to W 1 ’ 00 initial data 

In this step we consider a solution u of (2.8) granted by Theorem [2] and show 
that the requirement on H 2 smoothness of the initial data can be weakened to 


IT 1 ’ 00 . To this end we pass to the limit in (2.8) with an approximating sequence 
of smooth initial data. 

The following result establishes a bound on ||w||l°°(/) independent of the 
U 2 -norm of initial data (unlike in Lemma |4j) and provides also an estimate for 

IMIw 1 .°°(i)* 


Lemma 5. Let u be a solution of ( |2.8| ) (with initial value uq) on the interval 
[0,T] satisfying ||u(t)||z~(/) < Sq for all t <T. Then 


IK*)||z,°°(J) < ||wo||z°°(/) + Rt 


(2.38) 


where R > 0 is a constant independent of uq. Furthermore, the following in¬ 
equality holds 

||w<t||l“(/) < ai(t), (2.39) 

where ai(t) is the solution of 
da i 


dt 


— Pioi + QiOi + Ri, a(0) — || (i*o)o- II z,°° (z) 


(2.40) 


(continued by +oo after the blow up time) and Pi, Q\, R± are positive constants 
independent of Uq . 


Proof. Both bounds (2.38) and (2.381 are established by using the maximum 
principle. 

Consider g(a,t) = g(t) := ||uo||l“(/) + Rt, where R > 0 is to be specified. 
Since $ is bounded it holds that g t — < h(g t ) > R — sup>I>. Assuming that 
u — g attains its maximum on (0, t] x /, say at (oo,£o)> we have u a (cro,to) = 0, 
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u a<J {ao,t 0 ) < 0, g t (cr 0 ,to ) < u t {a 0 ,t 0 ) and S/(l - uk 0 ) = 1. Using ( |278[ ) and 
monotonicity of the function V — 4>(V) we also get for t = to, cr = cr 0 


„ * ,, > S /I — UKq \ 

R- sup$ < g t - <&{g t ) <u t - --$(--- u t ) 

1 — UK o V b / 


< 


«0 


Ut 

{l' i (^p Lu ’) Sd ° + 2,r )' 


1 — zi/«o L[u\ 


The last term in this inequality is uniformly bounded for all functions u satisfy¬ 
ing |u| < So, and therefore u — g cannot attain its maximum on (0, i] x I when 
R is sufficiently large, hence u(a, t) < ||uo||/,«>(/) +Rt. Similarly one proves that 
u(a,t) > -| KIU°° (/) - Rt- 


To prove (2.39) we write down the equation obtained by differentiating (2.8) 
$'(U) + 2 


in cr (cf. fl2.26D), 


'Mat 




(1 - $'(U))5 4 
S 


u a u a 


= Au n 


+ B+(- -) (d>(U) 

V 1 — UK 0 / cr V 


27T 


1 


L[u] L[u] 


<$>{V)Sda J, 

(2.41) 

where we recall that A(a, V, u,p ) and B(a, V, u,p) are bounded continuous func¬ 
tions on I x I x [— So, <5o] x K. Consider the function u (T (a 1 1 ) — ai(t), with a±(t) 
satisfying ai(0) = ||«<t( 0)||£°°(/). If this function attains its maximum over 
I x (0 ,t] at a point (cr 0 ,to) with t 0 > 0, we have at this point 


dai 

dt 


<u at <B 


S 

1 — UKq ) a ' 
< Cl + + CoU 0 




2 tt 

L [u\ 


1 

L[u] 


<5>(V)Sda 


(2.42) 


where C\, C 2 , C 3 are some positive constants independent of u. We see now 
that for Pi > C 2 , Q 1 > C 3 and Ri > C\, either u a (<T 0 ,to) < ai(io) or 
inequality (2.42) contradicts (2.40). This yields that for all cr £ / and all 
t < sup{r;ai(r) is finite}, 

u a (a, t) — ai(t) < max{0, max u a (a, 0) — ai(0)} = 0. 

o-e/ 

The lower bound for u a (cr,t) is proved similarly. □ 

Lemma [ 5 ] shows that for some T± > 0 the solution u of (2.8) satisfies 

||u||l°°(j) < and ||itcr||z,°°(/) < Mi, (2.43) 

when 0 < t < min{T, Ti}, where T is the maximal time of existence of u. 
Moreover, X) and constant Mi in (2.43) depend only on the choice of con¬ 
stants 0 < a < 1 and M > 0 in the bounds for the initial data |uo| < a and 
||(wo)<t||l”(/) < M. We prove next that actually T > T\. 
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Lemma 6 . Let u solve (2.8) and let (2.43) hold on 0 < t < T 2 , where 0 < X 2 < 
Ti. Then for any 6 > 0 


\ucia 1 , t) — u a {a", t )| < C&\a' — <j"\ when 6 < t < T 2 , 


(2.44) 


where 0 < d < 1 and C ,5 depend only on S and constant Mi in (2.43). 


Proof. Recall from Lemma[l]that 4' denotes the inverse function of V — 4>(R), 
and that 

so that the solution V of the equation V = 4>(R)+k—A is given by 4 , (k—A). This 
allows us to write the equation (2.41) for v := u a in the form of a quasilinear 
parabolic equation, 


v t = (a{v a ,a,t)) a - 
where 


\{t)i 


a(p, <r, t ) = 


S(u)( 1 — UKq) 

1 — UKq 


V a - X(t)(uK 0 ) c 


S(u)(l — UKq) 2 ’ 

($(4 '[R(p,a,t) - A(t)]) +R(p,a,t)j 


(2.45) 


with 


(p, <7, t) = 3 ^, ((1 - UK 0 )p + 2 K 0 ul + ( K 0 )crUcrU + K 0 (l - UK 0 ) 2 ^j . 


S 3 (u) 

Note that setting V := 4— X(t)] we have 


S(u) ( $'(V) \ dR. 

. -UKn(l-$'m + X ) &p( P,CT,t ^ ~ 


da 

dp 1 — un 0 V1 — 4>'(R) ' "J dp' 


1 


S 2 {u){ 1-$'(V))' 


It follows that if (2.43) holds then the quasilinear divergence form parabolic 
equation (2.45) satisfies all necessary conditions to apply classical results on 
Holder continuity of bounded solutions, see e.g. pt3] [Chapter V, Theorem 1.1]; 
this completes the proof. □ 

The property of Holder continuity established in Lemma [ 6 ] allows us to prove 
the following important result 


Lemma 7. Let u solve (2.8) and assume that (2.43) holds then for any r > 
5 > 0 


\ u crai.t)\\ 2 L2{I) < (||u ( 7 (T (r)||| 2 (J) + l)e P2(t T) when t> 


(2.46) 


where P 2 depends only on S and constant C in (2.43) (the latter constant depend 
not on the PI 2 norm of uq but on its norm in the space W 1,00 (I)). 
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Proof. Introduce smooth cutoff functions </>„(cr) satisfying 


</> n (cr) = 1 on [ 1 /n, 2 /n] 

0 < 4> n (&) < 1 on [0,1/n] U [2/n,3/n] and \(<j> n )a\ <Cn. (2.47) 
4 > n { cT ) = 0 otherwise, 


Consider t> 5. Multiply (2.41) by — (</> 2 u <jo) a and integrate in a, integrat¬ 
ing by parts in the first term. We obtain 

\jt Jj ^ nU ^ d(J + ^ Jj < t ) l U ™cr d °' < C J <t> 2 n{ula + W<J<r\ )\u aaa \do 

~\~ Cn J" (f n ( \ii(7(7(7 1 \ 'U‘CT(7 1 d - l^crcr | d - | ) do. 

(2.48) 

where 7 ' > 0 and C are independent of u and n. Applying the Cauchy-Schwarz 
and Young’s inequalities to various terms in the right hand side of ( 2.48| ) leads 
to 

\jt J <t>l u lcr da +^y J 4 > 2 n ul aa do < C'l J 4 > 2 n u A aa do+C 2 n 2 J (\u aa \ 2 +l)do, 

(2.49) 

where C'l, C 2 are independent of u and n and s(4> n ) denotes the support of 
4> n - Next we apply the following interpolation type inequality (see, e.g., EJ> 
Chapter II, Lemma 5.4) to the first term in the right hand side of (2.49): 


J 4 >n u ta d<J < C 3 ( sup{| u a (o') - u a {o ")|; o', o" G s( 0 „)})" 


X / <f>n U lcra d ^ + 


's(^n) 


! |(</>n)<7| 2 d<7). 


(2.50) 


Now we use (2.44) in (2.50) to bound sup-du^er') — u a (o ")|; o', 0 " G s(^>„)} 
by Cs^/ri)® and choose n so large that C 2 C 3 C 2 (l/n) 2,? < 1 / 47 ', then ( 2.49| ) 
becomes 


1 d r 

2 dt J j 


\ul„do + 


47 1 


,/ / t n 000 


u 2 do < C A n 2 


C(0„) 


4- l)do. (2.51) 


It is clear that we can replace </> n (cr) in (2.51) by its translations <p n (o + k/n), 
k G Z, then taking the sum of obtained inequalities we derive 

\j t J 4-^7 J $n U lav d ° < C 5 ( J \u aa \ 2 do + l) , (2.52) 

where C 5 is independent of u, and <p n = 0 2 (cr 4 -fc/n). Note that 1 <<(>„< 3, 

therefore applying Gronwall’s inequality to (2.52|) we obtain (2.46). □ 
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Corollary 1. Assume that a solution u of (2.8) exists on [0, T] for some T > 0, 
and (2.43) holds. Then, given an arbitrary positive 5 <T, we have 


l^fTererl dcdt ^ C$ || U acT (t) || jj2 , 


(2.53) 


for all 5 < t < T, where Cs depends only on 5. 

Proof. The bound follows by integrating (2.52) in time and using (|2.46[). □ 


Using Lemma [7] and Lemma [5j taking into account also Remark [7J we see 
that solutions in Theorem [2] exist on a common interval T = T(a,M), provided 
that ||u 0 || < a < 1 and || (w 0 ) < M. 

In the following Lemma we establish an integral bound for - 


Lemma 8. Let u be a solution (2.8) on [0, T] satisfying (2.43). Then 


l v<?\\L 2 (Ix[0,T}) 


<C, 


(2.54) 


where C depends only on the constants in (2.43). 


Proof. To obtain (2.54) one multiplies (2.8) by u aa integrates in er, integrating 


by parts in the first term. Then applying the Cauchy-Schwarz and Young’s 


inequalities and integrating in t, one derives (2.54), details are left to the reader. 

□ 

Various estimates obtained in Lemmas mm and Corollary [T] as well as 
Lemma[8]make it possible to pass to general initial data uq G WdJer° (/). Indeed, 
assume a := ||uo||z,°°(/) < 1 and let M := ||(ito)o-IIz,°°(z)- Construct a sequence 
Uq —^ u o converging weak star in Wp£?(I) as k —► oo, where u k G Hp er (I). 
This can be done in a standard way by taking convolutions with modifiers so 
that ||u§|| loo(j) < a and ||(wo)^IU°°(/) — M. Let u k (a,t) be solutions of ( |2.8| ) 
corresponding to the initial data u §(ct). We know that all these solutions ex¬ 
ist on a common time interval [0, T\ and that we can choose T > 0 such that 
( 2.43| ) holds with a constant C independent of k. By Lemma [8] the sequence 
u K (a, t) is bounded in L 2 (0,T; H 2 (I)). Therefore, up to a subsequence, u k (a,t) 
weakly converges to some function u(a,t ) G L 2 (0, T; H 2 (I)). Using (2.8) we 
conclude that u k (a,t) converge to u t (cr,t) weakly in L 2 (0, T; L 2 (J)). It follows, 
in particular, that ufc r, 0) = Uq(<t). Next, let 6 > 0 be sufficiently small. It fol¬ 
lows from Lemma|8l that there exists r G [5,25] such that ||u£ 0 .(t)|| < C/5. 
Then by Lemma [7^ and Corollary 0 norms of u k in L°°(25, T; H 2 er (/)) and 
L 2 (25, T; Hp er (I)) are uniformly bounded. Thus u k (a,t ) converge to u(a,t) 
strongly in L 2 (25, T; if 2 r (I)). This in turn implies that u k (a,t) converge to 


per \ 

u t {a,t) strongly in L 2 (0, T; L 2 (/)). Therefore the function u(a,t) solves (2.8) 
on [25, T]. Since 5 > 0 can be chosen arbitrarily small, Theorem[l]is completely 
proved. 
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3. Non-existence of traveling wave solutions 


The following result proves that smooth ( H 2 ) non-trivial traveling wave 
solutions of 11.41) do not exist. The idea of the proof is to write equations of 


motion of the front and back parts of a curve, which is supposed to be a traveling 
wave solution, using a Cartesian parametrization. Next we show that it is not 
possible to form a closed, H 2 -smooth curve from these two parts. We note that 
every traveling wave curve is always H 2 smooth since it is the same profile for 
all times up to translations. 

Theorem 3. Let $ satisfy conditions of Theorem^ IfT(a,t) is a family of 
closed curves which are a traveling wave solution of ( |1.4[ ), that is T(<r, t) = 
r(tr, 0) + vt, then v = 0 and T(cr, 0) is a circle. 


Proof. It is clear that if v = 0 then a circle is the unique solution of (1.4). 


Let T(cr,t) be a traveling wave solution of (1.4) with non-zero velocity v. By 
Theorem [T| T(-,f) is smooth ( H 2 ) for all t > 0. By rotation and translation, we 
may assume without loss of generality that v x = 0, v y = c with c > 0 and that 
T(cr, t) is contained in the upper half plane for all t > 0. Let T(cro, 0) be a point 
of r(cr, 0) which is closest to the x-axis. Without loss of generality we assume 
that r(cr o ,0) = 0. Locally, we can represent T(cr, t) as a graph over the x-axis, 
y = y{x ) + ct. Observe that the normal velocity is given by 


V = - 


v 71 + ( v '( x ) 2 ) 


and the curvature k is expressed as 


c(x) = 


y"(x) 


(:i + (y'(x )) 2 ) 3/2 ' 


(3.1) 


(3.2) 


Adopting the notation 

fx(z) ■■= 


\/o 


- $ 


it follows that y solves the equation 


where, by construction, 



(3.3) 

= fx(y') 

(3.4) 

^o 

II 

o 

(3.5) 


Observe that (3.4) is a second order equation for y which depends only on 


y'. Thus, we may equivalently study w := yf which solves 

w' = fx(w), w(0) = 0. 


(3.6) 


Note that (3.6) is uniquely solvable on its interval of existence by Lipschitz 
continuity of /£. Further, the definition of /_( guarantees that w has reflectional 
symmetry over the y-axis. 
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If (3.6 1 has a global solution, w, then it is immediate that y(x) := Jg w(s)ds 
cannot describe part of a closed curve. As such we restrict to the case where w 
has finite blow-up. 

Assume the solution wg of ( |3.6| ) has finite blow-up, wg( x) —> +00 as x —> x* B 
for some 0 < x* B < 00 . Then 9 b{x) := vj B (,s)d.s has a vertical tangent 
vector at ±x* B . To extend the solution beyond the point x* B , we consider wf, 
the solution of (3.6) with right hand side f^ c . As above we assume that wf 
has a finite blow-up at x* F > 0. Defining 9f{x) := f Q w F (s)ds, we have the 
following natural transformation, 


Vf{x) := -vf{x - ( x* B - x* F )) + 9b{x b ) + 9b{x* f ). (3.7) 


Note that gluing yp to ys forms an H 2 smooth curve at the point (x* B ,yB(x B )) 
if and only if w F (x) —► +00 as x —> x* F . We claim that this is the unique, 
smooth extension of ys at x* B . To that end, consider the rotated coordinate 
system (x,y) 1 —> (y,—x). In this frame, the traveling wave moves with velocity 
v x = c, v y = 0, and can be locally represented as the graph x = x(y) + ct, x(y) 
solving 

x"=9l{x') (3.8) 


with 


9l(z) ■■= 


—cz 


x/T 


- $ 


—cz 




+ A ) (l + z 2 ) 3/2 . 


(3.9) 


As before, g £ is Lipschitz and so solutions of (3.8) are unique, establishing the 
claim. 

To complete the proof, we prove that x F > x B , which guarantees that 
the graphs of 9b(x) and 9f{x) can not smoothly meet at — x* B . Due to the 
monotonicity of V — < b(I / ) we have that for any w, f%(w) > f^ c (w). Thus w' B > 
w' F for any fixed w. Since wb{ 0) = wf( 0), we deduce that wb{x) > w F {x) for 
all x > 0. It follows that x F > x B . Let X 2 £ (0, x^) and observe by continuity 
of wb that there exists X\ G (0, X 2 ) such that wb(x 1 ) = Consider the 

solution w of 

w' = fx(w) w{x 2 ) = wf(x 2 ). (3.10) 


Note that w(x) = wb(x — (x 2 — Xi)) and so the blow-up of w(x) occurs at 
x*b + X 2 — Xi. Since w(x) > w F (x) for all x £ (x 2 , x b +x 2 — X\) it follows that 
x F >x* B +x 2 — X\ > x* B . □ 
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4. Numerical simulation; comparison with volume preserving curva¬ 
ture motion 

The preceding Section [2] proves short time existence of curves propagating 
via (1.41, and Section [3] shows the nonexistence of traveling wave solution. In 
this section we will numerically solve (1.4) and for that purpose we will introduce 
a new splitting algorithm. Using this algorithm, we will be able to study the 
long time behavior of the cell motion by numerical experiments, in particular, 
we find that both non-linearity and asymmetry of the initial shape will result 
in a net motion of the cell. 

Specifically we numerically solve the equation Q written so that the de¬ 
pendence on p is explicit ( / 34> 0 (U) = $(V), see remark[l]): 


V(s, t ) = k(s, t ) + /34> 0 (U(s, t)) - 


\m\ 


(K(a',t)+P^ 0 (V(a',t)))ds' t (4.1) 


We propose an algorithm and use it to compare curves evolving by (4.1) with 
0 < (3 < P cr , with curves evolving by volume preserving curvature motion 
(/? = 0 ). 


For simplicity, we assume that 4> 0 (U) defined via (1.5) has a Gaussian shape 
4> 1 (U) := e“l v I" (see Figure [ij which agrees well with the actual graph of 4 >q. 
This significantly decreases computational time in numerical simulations. 


4-1- Algorithm to solve \4-A) 

In the case /3 = 0 (corresponding to volume preserving curvature motion), 
efficient techniques such as level-set methods pain] and diffusion generated 
motion methods nana can be used to accurately simulate the evolution of 


curves by (4.1). There is no straightforward way to implement these methods 


when P > 0 since V enters equation (4.1) implicitly. 


Moreover, due to the non-local term, a naive discrete approximation of (4.1) 
leads to a system of non-linear equations. Instead of using a non-linear root 
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V 


Figure 4: Graph of <t>o 


solver, we introduce a splitting scheme which resolves the two main computa¬ 
tional difficulties (non-linearity and volume preservation) separately. 

In particular, we decouple the system by solving the N local equations 

V i = K i +p$ 1 {V i )-C, (4.2) 


where C is a constant representing the volume preservation constraint which 
must be determined. For /3 < /3 cr , (4.21 can be solved using an iterative method 
which (experimentally) converges quickly. The volume constraint can be en¬ 
forced by properly changing the value of C . 

We recall the following standard notations. Let Pi = ( Xi,yi ), i = 1,..., N 
be a discretization of a curve. Then h := 1/N is the grid spacing and 


Dpi := 


Pi +1 — Pi -1 
2 h 


and D 2 pi := 


Pi+i - 2pi + pi_i 

h 2 


(h = 1/N) 


(4.3) 


are the second-order centered approximations of the first and second derivatives, 
respectively. Additionally, (a, &)-*- = (— b,a). 


We introduce the following algorithm for a numerical solution of (4.1). 


Algorithm 1. To solve (4.1) up to time T > 0 given the initial shape T(0). 
Step 1: (Initialization) Given a closed curve T(0), discretize it by N points p) = 
(•«)• 

Use the shoelace formula to calculate the area o/T(0): 

1 


A °=2 


x °iV°i+i 


r 0 ?/ 0 -’Vr 0 7/°-T°7/° 


(4.4) 


i—1 i=l 

Set time t := 0, time step At > 0, and the auxiliary parameter C := 0. 
Step 2: (time evolution) If t < T, calculate the curvature at each point, Ki using 
the formula 


det(Dpi D 2 p*) 

\mr 


(4.5) 
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where 

solve 


is the standard Euclidean norm. Use an iteration method to 


V? = Ki +£$!(■ V?)-C 


(4.6) 


to within a fixed tolerance £ > 0. 
Define the temporary curve 


where Vi = (Dp®) A /\\Dp®\\ is the inward pointing normal vector. 

Calculate the area of the temporary curve A temp using the shoelace formula 


(4.41 and compute the discrepancy 

AA := (A temp - A°) ■ (A 0 )- 1 , 

If | AA| is larger than a fixed tolerance e, adjust C i—>• C + A A and return 


to solve (4.6) with updated C. Otherwise define p\ +At \=p*f rnp and 


T(t + At) :={p* +At }, 

Let t = t + At, if t <T iterate Step 2; else, stop. 

In practice, we additionally reparametrize the curve by arc length after a 
fixed number of time steps in order to prevent regions with high density of 
points pi which could lead to potential numerical instability and blow-up. 


Remark 8. In Step 2, for /3 < 1, the right hand side of (4.6) is contractive 


and thus we may guarantee convergence of the iterative solver to the solution of 


(4.6). 


We implement the above algorithm in CH —V and visualize the data using 
Scilab. We choose the time step At and spatial discretization step h = — so 
that 

At 1 
h? - 2 ’ 

in order to ensure convergence. Further, we take an error tolerance 

e = .0001 


in Step 2 for both iteration in (4.6) and iteration in C. 


j.2. Convergence of numerical algorithm 

To validate results, we check convergence of the numerical scheme in the 
following way: taking a fixed initial curve F(0), we discretize it with varying 
numbers of points: N = 2 m for m = 5,..., 8, and fix a final time step sufficiently 
large so that T(kAt) reaches a steady state circle. Since there is no analytic 
solution to (1.4), there is no absolute measure of the error. Rather, we define 


the error between successive approximations N and 2 N. 

To this end, we calculate the center location, Cn = (Cj^, C'jy), of each steady 
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state circle as the arithmetic mean of the data points. Define the error between 
circles as 

errjv = ||Cjv — C 2 N\\p- 

Then, the convergence rate can be expressed as 

p := lim pn = lim log 2 — ——. 

N—>oo N—>oo err 2 at 

We record our results in Table]!] as well as in Figure[5] Note that p ~ 1, namely 
our numerical method is first-order. 


N 

err/v 

Pn 

32 

.2381754 

.9141615 

64 

.1263883 

.9776830 

128 

.0641793 

1.0019052 

256 

.0320473 

- 


Table 1: Table for convergence of numerical simulations 



Figure 5: Convergence of numerical simulation for decreasing mesh size h = 1/N (zoomed in) 


4.3. Numerical experiments for the long time behavior of cell motion 

We next present two numerical observations for the subcritical case /3 < (3 cr . 
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1 . If a curve globally exists in time, then it tends to a circle. 

That is, in the subcritical /3 regime, curvature motion dominates non-linearity 


due to d>(V). This is natural, since for small (3 the equation (4.1) can be viewed 


as a perturbation of volume preserving curvature motion, and it has been proved 
(under certain hypotheses) that curves evolving via volume preserving curvature 
motion converge to circles (8]. 


In contrast, the second observation distinguishes the evolution of (4.1l from 
volume preserving curvature motion. 

2. There exist curves whose centers of mass exhibit net motion on a finite 
time interval (transient motion). 

The key issue in cell motility is (persistent) net motion of the cell. Although 
Theorem [3] implies that no non-trivial traveling wave solution of (1.4) exists, 
observation 2 implies that curves propagating via (4.1) may still experience a 


transient net motion compared to the evolution of curves propagating via vol¬ 
ume preserving curvature motion. We investigate this transient motion quanti¬ 
tatively with respect to the non-linear parameter /3 and the initial geometry of 
the curve. 


4.3.1. Quantitative investigation of observation 2 

Given an initial curve T(0) discretized into N points p ° and given 0 < /3 < 
/3 cr , we let T^fcAf) := {p^} be the curve at time fcAi, propagating by (4.1). 
In particular, To(kAt) corresponds to the evolution of the curve by volume 
preserving curvature motion. 

Our prototypical initial curve T(0) is parametrized by four ellipses and is 
sketched in Figure [bj 
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(4cos(0), 3sin(0)) for — ^ < 9 < ^ (4.7) 

0 7T Q 77 - 

(2cos(0), - sin(0) + -) for - < 9 < — (4.8) 

(C cos(6>), ^ sin(6>)) for - | < 9 < | (4.9) 

Q 7T Q’7T 

(2cos(0), - sin(0) - -) for - < 9 < (4.10) 

The parameter Q determines the depth of the non-convex well and is used as 
our measure of initial asymmetry of the curve. 



Figure 7: Overall drift of curve with = 1 (blue, solid) compared to /3 
starting from initial data (4.7)-(4.10) with £ = 2 


0 (green, dashed), 


To study the effect of /? and asymmetry on the overall motion of the curve, 
we measure the total transient motion by the following notion of the drift ofTp. 
First fix T(0). Note that observation 1 implies that for sufficiently large time 
kAt, Tp(kAt) and To(fcAt) will both be steady state circles. Define the drift of 
Tp to be the distance between the centers of these two circles. 

Remark 9. Note that this definition is used in order to account for numerical 
errors which may accumulate over time. Numerical drift of the center of mass 
of the curve caused by errors/approximations is offset by “calibrating” to the 
/3 = 0 case. 


We consider the following two numerical tests: 

1. Dependence of drift on (3: Starting with the initial profile (4.7)-(4.10) with 
C = 1, compute the drift of Tp for various values of /3. 
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2. Dependence of drift on initial asymmetry: Starting with the initial curve 
with /3 = 1, compute the drift of Yp for various values of f. 

Taking T = 20 is sufficient for simulations to reach circular steady state. We 
observe that drift increases with respect to f and increases linearly with respect 
to /3. These data are recorded in Figure [8] 


(4.7l-(4.10 



Figure 8: Dependence of drift on the parameters /3 and £ 
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